Elliptic Flow from a Transversally Thermalized Fireball 
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The agreement of elliptic flow data at RHIC at central rapidity with the hydrodynamic model 
has led to the conclusion of very rapid thermalization. This conclusion is based on the intuitive 
argument that hydrodynamics, which assumes instantaneous local thermalization, produces the 
largest possible elliptic flow values and that the data seem to saturate this limit. We here investigate 
the question whether incompletely thermalized viscous systems may actually produce more elliptic 
flow than ideal hydrodynamics. Motivated by the extremely fast primordial longitudinal expansion 
of the reaction zone, we investigate a toy model which exhibits thermalization only in the transverse 
directions but undergoes collisionless free-streaming expansion in the longitudinal direction. For 
collisions at RHIC energies, elliptic flow results from the model are compared with those from 
hydrodynamics. With the final particle yield and k± -distribution fixed, the transversally thermalized 
model is shown not to be able to produce the measured amount of elliptic flow. This investigation 
provides further support for very rapid local kinetic equilibration at RHIC. It also yields interesting 
novel results for the elliptic flow of massless particles such as direct photons. 

PACS numbers: 25.75.-q, 24.10.Nz, 25.75.Ld, 24.85. +p 



I. INTRODUCTION 



Recently the STAR collaboration showed jl], |2|, || that 
the k± and centrality dependence of the measured elliptic 
flow coefficient V2 of pions and protons at the Relativis- 
ts Heavy Ion Collider (RHIC) largely agreed with those 
generated from the hydrodynamic model simulation of 
heavy ion collisions |, § | §. This was subsequently 
confirmed by data from the PHENIX || and PHOBOS 
H collaborations. V2 is the second coefficient of an az- 
imuthal Fourier expansion of the transverse momentum 
spectrum around the beam axis ]lO| and, for collisions be- 
tween identical nuclei, is the lowest non-zero anisotropic 
flow coefficient at midrapidity. Given that the main pre- 
requisite of the hydrodynamic model is complete local 
thermal equilibrium, which requires very intense rescat- 
tering among the matter constituents, it was argued that 
hydrodynamics should give the largest possible elliptic 
flow, and the observation that the data saturate the hy- 
drodynamic limit was taken as evidence that thermal- 
ization must be very fast at RHIC for central and semi- 
central collisions || ^, . 

Elliptic flow requires reinteractions within the pro- 
duced matter as a mechanism for transferring the initial 
spatial deformation of the reaction zone in non-central 
collisions onto momentum space. It it thus plausible to 
expect that the largest elliptic flow signal is produced in 
the hydrodynamic limit, i.e. in the limit of infinite rescat- 
tering rates (|, ^, |ll| ; however, a proof for this hypothesis 
has not yet been found. In hydrodynamics, collective flow 
is generated by pressure gradients, and flow anisotropics 
such as elliptic flow require anisotropic pressure gradi- 
ents in the plane transverse to the beam direction. The 
longitudinal pressure plays a less visible role: the work 
done by the longitudinal pressure reduces the transverse 
energy dEx/dy and thereby, for given initial conditions 



for the energy density deposited in the reaction zone, the 
amount of transverse flow at freeze-out. If there were no 
longitudinal pressure, more of the initially deposited en- 
ergy would go into the transverse directions, leading to 
flatter transverse momentum spectra. It is, at least in 
principle, conceivable that this could also lead, at fixed 
radial flow, to larger anisotropies of the transverse col- 
lective flow. If this were indeed the case and, at fixed 
slope of the angle-averaged spectra, larger values of V2 
could be generated in this way than within the usual 
hydrodynamic, the early thermalization argument would 
break down since the RHIC data [jl] |2|, [| would no longer 
saturate the theoretical prediction from such a model. 

To explore this hypothetical possibility, we study in 
this paper a toy model which assumes that, at high colli- 
sion energies, during the earliest collision stages the longi- 
tudinal momenta of the produced particles do not ther- 
malize, but that the strong initial longitudinal expan- 
sion is instead dominated by collisionless free-streaming. 
The initial transverse momenta are much smaller and 
assumed to thermalize quickly. This results in a sys- 
tem with local transverse, but zero longitudinal pressure. 
Due to the masslessness of the partons created in the re- 
action zone the trace of the energy momentum tensor 
vanishes; in the absence of longitudinal pressure Py , the 
transverse pressure Pj_ must thus be related to the initial 
energy density e = T 00 by e — 2Pj_ (instead of the usual 
e = 3P) . This results in a much stiffer equation of state 
for the transverse dynamics (see Sec. II B), giving more 
flow and possibly stronger flow anisotropies. Of course, 
if the transverse momenta thermalize, the longitudinal 
momenta should do so eventually, too. We assume in our 
model that when this happens the flow anisotropies have 
already almost reached their asymptotic values Jfjj. 

Similar to the hydrodynamic simulations ||, |[ [r|, 
we assume that the longitudinal free-streaming dynam- 
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ics is boost invariant H] and correspondingly concen- 
trate on the central rapidity region. Starting from a ki- 
netic description of the longitudinally free-streaming but 
transversally thermalized gluonic system (which we call 
"transversally thermalized model", TTHM), we derive 
a set of macroscopic evolution equations (TTHM equa- 
tions) which we solve with appropriate initial conditions. 
The impact parameter dependence of the initial condi- 
tions is handled in the same way as in the hydrodynamic 
simulations. When comparing our TTHM solutions to 
those from the ideal hydrodynamic model (HDM) we re- 
tune the initial conditions in such a way that for central 
collisions (6 = fm) roughly the same multiplicities and 
spectral slopes at midrapidity are obtained. We then 
compare the momentum-space anisotropy V2 in the two 
approaches, as a function of transverse momentum k± 
and impact parameter b. 



II. CONSTRUCTING A TRANSVERSALLY 
THERMALIZED SYSTEM 

A. The Phase-Space Distribution Function 

We consider a gluon dominated system at a short time 
after nuclear contact. We assume that at that time the 
transverse gluon momenta are already thermalized but 
that the system is still free-streaming along the beam 
direction. The system is required to possess longitudi- 
nal boost invariance, reflecting a boost invariant primary 
particle production mechanism For easiest imple- 

mentation we assume that the collision energy is so high 
that the colliding nuclei are Lorentz contracted to two 
infinitesimally thin sheets in the z-direction and that all 
produced particles point back to z = t = 0. Their longi- 
tudinal momenta and coordinates thus satisfy z = u z t or 
koz = k z t where u z — k z /ka. This automatically identi- 
fies Q the rapidity 



y 2 \ko-k z J 



(2.1) 



'k Q + k 

of the produced particles with their spacetime rapidity 



' 2 \t-z 



(2.2) 



We thus make the ansatz 



f(k, x, t) = S(k z - k z t) h(k, x, t) . (2.3) 

Longitudinal free-streaming is implemented by the con- 
dition 



Since the 5-function factor solves this equation trivially, 
this implies a constraint for h. The latter is most conve- 
niently expressed by using an alternate coordinate sys- 
tem (we hereafter make explicit use of the masslessness 



of the p roduced gluon degrees of freedom, by setting 
m± = y/ m 2 + k\ — > k±): 

t = t cosh 77 , z = rsmhi] , (2.5) 

fco = fc^coshy, k z = k± sinh y . (2.6) 



In these coordinates the ansatz Eq. (2.3) becomes 
5(r)- y) 



f(k,x,t) 



k±r cosh(?7 — y) 



h(k±,x±,r)-y,T) (2.7) 



where the dependence of h on only k± and the difference 
rj — y results from the requirement of longitudinal boost 
invariance. We implement transverse thermalization by 
the ansatz 



h(k±,x±,r)-y,r) 



Tpk± 
7-L 



cosh(rj—y) 



1 



jfc-u/T _ I ' 



(2.8) 



where u M =7^(cosh?y, v±. sinh 77) is the flow vector with 
Bjorken longitudinal flow velocity v z = tanh?7 [|14| and 



transverse flow velocity v± = v±(x±, t) (7J 1 = \JT—v\). 
T = T(xj_, t) is the temperature characterizing the ther- 
malized transverse momentum spectrum. In the local 
rest frame T is a function of the rest frame coordinates 
x* ± and t* — t* which after the standard Jl5| succes- 
sive transverse and longitudinal boosts to the laboratory 
frame, with vj_ and v z = tanh rj respectively, are related 
to the laboratory variables by 



t* = 



<<X|| 



7_i_(t cosh 7] — z sinh r/ — v± ■ x±) 
j±(t — v± ■ x±) = u ■ x (2.9) 
7_l_ {x± 11 — v± (t cosh rj — z sinh rj) ) 

7j_(a;^H — v±t) = —m ■ x (2-10) 

x ±± . (2.11) 



The additional || and _L subscripts denote the compo- 
nents of x± parallel and perpendicular tovj, respec- 
tively, and the vector m M is defined in Eq. ( 2.16 ) below. 
The sequence of boosts described here agrees with the 
standard choice fa] and gives the simplest structure for 
the set of 4-vectors needed for de comp osition of the en- 
ergy momentum tensor (see Sec. [IB). As one can see 
from Eqs. (2.£)-(2.11), the dependence of T on x* ± and 
t* in the local rest frame translates into a dependence on 
x± and r in the laboratory frame. 

By virtue of the factor Sjrj —y), the argument of the 
Bose distribution in Eq. (2.8) reduces in the local rest 
frame to k-u — k±, showing explicitly that the sy stem 
is only transversally thermalized. Combining Eqs. (2.7) 
and ( |2.8| ) and exploiting the 5-function, the distribution 
function simplifies to 



f(k,x,t) 
h(k±,x±,t) 



S(r) -y)h(k±,x±,t) 



(2.12) 
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In the new coordinates, / satisfies the relativistic trans- 
port equation 



B. Energy-Momentum Tensor and Conservation 
Laws 



(k± coshf?]— y) — — sinh(r] 

V or t 



-y, 



d_ 

dr/ 



+ k±-Vx)f{k,x,t) = C{k,x,t), (2.13) 

where the collision term C is responsible for keeping the 
transverse momenta thermalized while not changing any 
of the longitudinal momenta. Due to the implemented 
free-streaming properties we have 



d 1 d\ f t 

cosh(?7-y)— - sinh^-y)-— I I — 8(r) - y) ) = 



such that Eq. ( 2.13| ) reduces to the following kinetic equa- 
tion for the transverse distribution function: 



fcL^ + fcx-Vx S Mk_.x ± .r) = C(k.x.l). (2.1.4! 



The space-time evolution of h(k±,x±,r) is entirely due 
to the collisions and collective transverse expansion, but 
decoupled from the boost invariant longitudinal expan- 
sion. The collisions must be sufficiently dominant over 
the transverse expansion to maintain the equilibrium 
form of h(k±, x±, r) given above. 

We are aiming at a macroscopic description on the ba- 
sis of the differential conservation laws for energy and 
momentum, similar to hydrodynamics. This requires 
the construction of T^ v from the distribution functions. 
Due to the different microscopic physics in the longitudi- 
nal and transverse directions implemented in our model, 
there are now more vectors required to construct a com- 
plete set of tensors with respect to which T^ v should be 
decomposed. The set of four- vector fields that we need 
are 

u tl (x± 1 ri,T) = 7x(cosh?7, vx, sinhry), 
n ,J, (x±,?],T) = (sinhry, 0, 0, cosh^), (2.15) 
m^ccx, r], t) — 7x(v_i_ cosh 77, v±, v±_ sinh?7). 

They arise from the following local rest frame vectors by 
successive boosts with v± and v z as described above: 

u"(asj.,»7,T) = (1,0,0,0), 
n»(x±,t},T) = (0,0,0,1), 
m^(x ± ,r],T) = (0,t>j.,0). (2.16) 

One is timelike, u 2 = 1, while the other two are spacelike, 
n =m =— 1. They are mutually orthogonal: u-n — 
u-m = n-m = 0. With these vectors, the distribution can 
be expressed in a form which is explicitly boost invariant 
under longitudinal boosts: 

mx,t) = ^ s ^- n ^^ T -_ v - k - m K (2.i7) 



Based on the structure ( 2.17 ) for the phase-space dis- 
tribution function we want to derive a macroscopic set 
of dynamical equations based on the conservation laws 
for energy and momentum. In terms of the distribution 
function /, the energy- momentum tensor is given by 



T" v (x,t) 



d 3 k k^k u 
(2tt) 3 k° 



f(k,x,t) 



(2.18) 



where v g = 2 • 8 = 16 is the degeneracy factor for color and 
helicity of the gluons. From the previous subsection we 
know that f(k,x,t) depends on the three 4-vectors tt M , 

and m M . The most general form of the tensor 
thus reads 

T iw = Au^u" + Bn^n v + Cm^m v + D^r^+n^u") 
+F(u fl m' J +m, fl u 1 ') + G(n^m" '+m tl n v ) + Hg^ (2.19) 

Since the gluons are massless, the trace of T^ v vanishes: 



rp/j, 



A - B - C + 4ff 
In the local rest frame we have 







T *o = A 



H = e, 



xx yy 
T* = B-H = P Z , 



(2.20) 



(2.21) 



where e is the energy density and Pj_ and P z are the 
transverse and longitudinal pressure, respectively. In 
a ge neral frame, all the Lorentz scalar coefficients in 
Eq. ( 2.1£ ) can be found by contracting with all pair- 
wise combinations of the 4-vectors i// 1 , n M and m^. We 
will denote these contractions by (uTn) = u^T^n^ etc. 
The details of this calculation are given in Appendix |A|) 
where we find 



(uTu) 
(uTm) 

(mTm) 

(uTn) 
(nTn) 
(nTm) 



A + H 
-F 

C-H 

-D 

B-H 

G 



T 

V, 



9 ~6fT 



TO 7T T 
— Vi 
T 



TO 
T 



■ 

■ 



t - 120 
-v„ 



120 



1 

2 £ 
P 2 



From this and Eq. ( 2.21 ) we see that 
B = H = -Pl , A = e + 
and the tracelessness condition gives 



Pi 



£ + Pl 



Pl 



1 



3Px 



1 



Pi 



(2.22) 
(2.23) 

(2.24) 

(2.25) 
(2.26) 
(2.27) 

(2.28) 
(2.29) 
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Therefore the equation of state is 
e = 2Pi . 



(2.30) 



Combining this with the expression for (mTm) one finds 
C = 0. The expression for the energy momentum tensor 
is then 

= (e + PiJuV - Px_n»n v - P±g fiv 

+P ± v ± {u^m y + m^u") . (2.31) 

For completeness and later use we also give the expres- 
sion for the gluon density in the local comoving frame: 



. r C(3)T 3 
n g = u-j = -v g —^, 

where j^{x) is the gluon number current: 
d 3 k k» 



f{x,t) = v g 



(2tt) 3 k° 



f(k,x,t) . 



(2.32) 



(2.33) 



We now proceed to derive the macroscopic evolution 
equations. We first rewrite the conservation laws 



dvT w = o 

in (r), t) coordinates, 



(2.34) 



cosh 77 



V di 



drj 



sllm?? ( __ + ___ j 

(2.35) 



Vj_T v = , 



where the sum over i goes over the two transverse direc- 
tions. When writing these out explicitly in terms of e, 
Pi, and the longitudinal and transverse flow velocities, 
using Eqs. (2.31) and (2.16), one finds that all explicit de- 
pendence on 77 disappears. Boost-invariant initial condi- 
tions (i.e. ^-independent initial expressions for e and P± 
and the structure ( 2.16 ) for the vectors u M , n M , and m^) 
are thus preserved by the macroscopic evolution equa- 
tions. 

For the further analysis we can thus concentrate on 
the dynamics in the central transverse plane at z = rj = 0. 
There we have 



u" = 7x(l,«x,0), 
n» = (0,0,0,1), 
TO M = -f±{v ± ,v ± ,0) . 



T x(0, 0,0,1), (2.36) 



drj 

^- = (1,0,0,0), (2.37) 
' =7±(0,wl,«lH2.38) 



drj 



Using the ry-independence of e and P± we find 



drj 



= 7 i(e + P x )(5 0fl + v l L 8^) - P ± 5 0fi 



+~/lP± (2vl6»° + (1 + vDv^S^ ' 
T oo 6 o^ + T oi § iii = T o M ; (2.39) 



and the equations of motion (2.35) at z = rj - 



dr 

Q T 0t 

~lh~ 



+ V 3 ,T^ = 0, 



V 3 T Jt = 0. 



= reduce to 
(2.40) 
(2.41) 



Using the equation of state (2.30) to eliminate e, these 



three equations completely determine the three unknown 
functions P± and v±. Since gluons possess no conserved 
charges, no further current conservation laws need to be 
considered. 



III. HYDRODYNAMICS 

For later comparison we here shortly review the paral- 
lel procedure for ideal hydrodynamics based on a locally 
fully thermalized phase-space distribution function 

In this case the energy-momentum tensor has the ideal 
fluid decomposition 

I"*" = (e + PjlAt" - Pg^ , (3.2) 

and for a longitudinally boos t-invariant system u M has 
the same form as in Eq. ( 2.16 ). The equation of state is 
in this case 

7T 2 T 4 

e = 3P = v„ 



30 ' 

and the gluon density in the local rest frame is 

C(3)T 3 



7T^ 



(3.3) 



(3.4) 



Note that the the energy per particle 
e 7T 4 

— = — -—Tk2.1T (3.5) 
n g 30C(3) V ; 

agrees exactly wit h th e corre sponding expression in the 
TTHM, see Eqs. (|2~22|) , (|2~32|) . 

The pressure is now locally isotropic. With longitudi- 
nally boost-invariant (77-ind cpcn dcnt) initial conditions 
for e and P the ansatz Eq. ( [2.16| ) for the flow velocity is 
preserved in time by the equations of motion d v T^ LV = 
H, and the latter become 



dr 
dT oi 



on 



T 



V J , T 3i 



P 



dr 



V J ± T ji = 



(3.6) 
(3.7) 



The main difference between these hydrodyna mical equa- 
t ions and our TTHM toy model equations ( 2.40] ) and 
(|2.41 ) is the pressure term on t he r ight hand side of 
Eq. ( 3~6| ) which is absent in Eq. ( 2.40p . It indicates the 
longitudinal work done by the isotropic pressure in hy- 
drodynamics; in the transversely thermalized (TTHM) 
model there is no longitudinal pressure which could do 
work in the longitudinal direction. 



IV. INITIAL CONDITIONS 

We consider Au+Au collisions at RHIC energies. The 
initial energy density distribution in the transverse plane 
as a function of impact parameter b is taken to be pro- 
portional to the density of binary nucleon collisions in the 
transverse plane, calculated from the normalized nuclear 
thickness function: 

T A {x x ) = J p A (x ± ,z)dz, J d 2 x ± T A {x ± ) = A. (4.1) 

For the nuclear number density we took the Wood-Saxon 
parametrization 



0.12 



Pa(x ± ) = 



;(r-Ho)/€ + 1 



(4.2) 



with R = 1.14 ^l 1 / 3 fm and £ = 0.54 fm. For collisions 
between nuclei of masses A and B the number of binary 
nucleon collisions per unit area in the transverse plane is 



dN(x ± ,b) 
d 2 x± 



a T A (x ± )T B (b- Xl _) 



(4.3) 



where Co = 40 mb is the total inelastic nucleon- nucleon 
cross section. We assume that the generated energy den- 
sity is proportional to this density. We performed two 
classes of calculations: in the first class we used the same 
maximum initial energy density as the authors of [[jj |7j 
(eo = eo = 23.0 GeV/fm 3 for central (6 = 0fm)) collisions 
at an initial time tq = 0.6 fm/c) in order to compare the 
HDM and TTHM results for similar initial conditions. In 
the second class of simulations we retuned the initial con- 
ditions in TTHM such that the same multiplicity density 
dN/dy and the same slope for the transverse momentum 
spectrum at midrapidity as in HDM is obtained. Both 
sets of results will be discussed in the next Section. 



V. ELLIPTIC FLOW IN A TRANSVERSALLY 
THERMALIZED SYSTEM 

A. Freeze-out prescription 

Just as the HDM, our TTHM model describes the 
time evolution of macroscopic thermodynamic quantities 
which must be converted to particle spectra before one 
can compare with experiments. To do so we use the 
well-known Cooper- Frye prescription pl| which gives the 
particle spectrum in terms of an integral over the phase- 
space distribution function along a so-called freeze-out 
hypersurface S(a;): 



E 



dN 



dN 



d 3 k dyk ± dk ± d(/) (2tt) 3 



k-d 3 T,(x) f(x,k) . 



, (5-1) 

d 3 £ M (x) is the four- vector integral measure normal to the 
hypersurface, and in the TTHM case f(x, k) is a transver- 
sally thermalized distribution function of the form (2.12), 
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FIG. 1: V2 as a function of impact parameter at two freeze- 
out energy densities gf . The solid lines are from our TTHM 
model while the dashed lines are from hydrodynamic simula- 
tions with identical initial conditions. 



with flow velocity u M evaluated along the freeze-o ut hy - 
persurface S and temperature T calculated via Eq. ( 2.22| ) 
from the energy density on S. Longitudinal boost- 
invariance dictates freeze-out along a surface of fixed lon- 
gitudinal proper time Tt(x±), and we can write 



S M (x) = (t{(x±_) cosily, x±, Tf(x±) sinhr/). (5.2) 
This gives 



k ■ d 3 E = ( k± cosh(j/-ry) - kj_ ■ Vj_T f J r f drj d z x± . (5.3) 

Of course, the ^-integration is trivi al du e to the factor 
5{rj—y) in the distribution function ( 2.12|) . 



B. TTHM vs. HDM for identical initial conditions 

In Fig. |l| we compare results for the elliptic flow from 
the TTHM and HDM models for identical initial con- 
ditions, taken from ^, Since the temperature pa- 
rameter T has a different meaning in our model than 
usual (only the transverse momenta are thermalized), we 
enforce "freeze-out" (i.e. we stop the dynamical evo- 
lution and calculate the spectra and momentum space 



TABLE I: Rapidity densities at y — and final times Tf for 
central collisions with identical initial conditions (see text) 
in the TTHM and HDM simulations. 





TTHM 


HDM 


E f (GeV/fm 3 ) 


dN/dy n (fm/c) 


dN/dy t { (fm/c) 


0.6 


780 6.95 


944 5.97 


0.8 


780 6.48 


940 5.36 



anisotropies) along a surface of constant energy density 
£f . We show results for two values of this parameter both 
of which lie in the quark-hadron transition region. Since 
hadronization certainly involves longitudinal momentum 
exchange, it will violate our assumption of longitudinal 
free-streaming, and we should definitely not follow the 
TTHM dynamics beyond the hadronization point. Of 
course, by truncating the dynamics at such high freeze- 
out energy densities we forfeit the possibility of compar- 
ing our results directly with experiment. However, a large 
fraction of the elliptic flow should already have developed 
at this point || 0, and we can still make a meaning- 
ful comparison between the hydrodynamic evolution of a 
fully thermalized quark-gluon plasma and that of a longi- 
tudinally free-streaming gluonic system with only trans- 
verse thermalization, by comparing the gluon spectra at 
this common "final" energy density £f. 

Figure [l] shows «2 as a function of the impact param- 
eter b. v%(b) = (cos(2</>))(6) is computed from the k±- 
integrated gluon spectrum at freeze-out as 



v 2 (b) 



J k±dk±d<p cos(20) 



dN 



dy k±dk±d(fi 



(b) 



J k±dk± 



dy k± dkj_d<j) 

(I tjilo 



(b) 



(5.4) 



The Figure shows that, for identical initial conditions, 
the TTHM dynamics generates a larger momentum-space 
asymmetry vi than hydrodynamics. We also see that v 2 
has not yet fully saturated at these values of ef, the reason 
for this can be found in Table | which shows that, with 
the very hard equation of state and initial conditions used 
here, the above values of £f are reached quite early, before 
the initial spatial deformation has been fully eliminated 
(see also Fig. 1 in Q ) . 

The larger values of v-i from TTHM, however, come 
also with much flatter transverse momentum spectra, 
as shown in Fig. |2[ Part of this flattening is due to a 
higher freeze -out temp erat ure in the TTHM. By com- 
paring Eqs. ( |2.22[ ) and ( |3.3| ) one sees at the same energy 
density £f the temperatures in the two models are related 
by T f , T = (2r f /To) 1/4 T f3 . With the r f values given m 
Table |, the TTHM freeze-out temperature is about a 
factor 2 higher than the HDM one. However, the inverse 
slopes of the spectra shown in Fig. ^ differ by more than 
this factor 2 (the difference is closer to a factor 2.5-2.6). 
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FIG. 2: Gluon k± distribution for central collisions (b = Ofm) 
and two different values for the freeze-out energy density, for 
identical initial conditions. The solid lines are for our model 
while the dashed lines are for hydrodynamics. 



This implies that TTHM also creates somewhat stronger 
radial flow than HDM which, in view of its harder equa- 
tion of state, is not unexpected. 

The initial conditions for the hydrodynamic model 
(in particular the initial energy density) were tuned Q 
to reproduce the measured final multiplicity densities 
dN c h/dy and spectral slopes in central Au+Au collisions 
at RHIC; clearly, with these same initial conditions, the 
transversally thermalized model TTHM will no longer be 
anywhere close to these data. The larger v 2 values from 
TTHM in Fig. [I] therefore cannot be compared with the 
RHIC data, and the comparison with HDM in Fig. [j] is 
misleading. For a meaningful comparison the initial con- 
ditions for TTHM should first be adjusted in such a way 
that at least the multiplicity densities and spectral slopes 
for central collisions arc similar in both models. This will 
be our next step. 

Before doing so, let us shortly comment on the ra- 
pidity densities dN/dy and collision durations listed in 
Table |[ For identical initial conditions, the TTHM gives 
considerably lower multiplicity densities. This is mostly 
due to the higher freeze-out temperature in the TTHM 
which implies that each gluon carries more energy (see 
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FIG. 3: The freeze-out time T{(x, y) as a function of the trans- 
verse coordinates x and y, respectively, for a freeze-out energy 
density S{ — 0.6 GeV/fm 3 . The solid (dashed) lines are for 
TTHM (HDM), respectively. The upper set of curves is for 
central collisions (6 = fm) where the freeze-out surface is az- 
imuthally symmetric. The lower set of curves corresponds to 
b = 7fm where we present cuts along the x and y axes, as indi- 
cated, to show the spatial azimuthal anisotropy at freeze-out. 
The impact parameter b points in x direction. 



Eq. (3J5)), so the same freeze-out energy density £f trans- 
lates into fewer gluons. This argument overestimates the 
difference between the models, though; the actual dif- 
ference is smaller since in the HDM the calculation of 
dN/dy involves an integration over space-time rapidity 
rj which can be interpreted as taking an average over 
fireballs with different flow rapidities rj ^ y and corre- 
spondingly reduced (redshifted) effective temperatures 
T eff = T f /cosh(?7-y). 

Longitudinal boost-invariance implies that the total 
transverse energy per unit rapidity, dEx/dy, is inde- 
pendent of y. Due to the absence of longitudinal pres- 
sure, no longitudinal work is done in the TTHM and 
dEr/dy is a constant of motion. This is different from 
the HDM where work done by the longitudinal pres- 
sure reduces dEr/dy with time. (Both statements were 
checked numerically.) Table | shows that, within numer- 
ical accuracy, the gluon multiplicity per unit rapidity, 
dN/dy, is also constant in time. While in ideal hydro- 
dynamics (HDM) this is a simple consequence of en- 
tropy conservation (for boost-invariant systems dS/dy 
is a constant of motion E3), and for massless particles 
the entropy per particle is independent of temperature), 
there is no such simple reason for this observation in 
the TTHM which is far from local equilibrium. De- 
tailed checks revealed that, within numerical accuracy, in 
the TTHM not only the transverse energy per particle, 
E T /N=(dE T /dy)/(dN/dy), but in fact the entire k±- 
spectrum of the gluons is completely time-independent! 



If the system were transversally homogeneous and ex- 
panded only longitudinally, this would not be surprising: 
in the a bsenc e of transverse gradients the conservation 
law Eq. (|2.40D for T 00 simplifies to 



d(rT 



00 \ 



0t 







(5.5) 



with T 00 = e cosh 2 i). In TTHM a consta nt pr oduct er 
implies a constant temperature (see Eq. ( 2.22| )) which, 
in the absence of transverse flow, then leads to a time- 
independent transverse momentum spectrum. In our 
case, however, the system expands in the transverse di- 
rection and cools. Our studies show that, in the longitu- 
dinal rest frame, the resulting loss in thermal transverse 
energy per particle is exactly compensated by the gain 
in transverse collective motion energy, in a way which 
exactly preserves the shape of the transverse momentum 
distribution. Although one should think that there must 
be a simple reason for this intriguing behaviour, we have 
not been able to find a simple analytical proof and can 
thus only present the numerical evidence. 

The longer freeze-out times in the TTHM simulations 
(see Table |) are a reflection of the lack of longitudinal 
pressure. In the HDM the pressure performs longitudinal 
work and thus causes a more rapid decrease of the energy 
density with time than in the TTHM. This is further 
elucidated in Fig. || where we show the freeze-out time 
Tf(x, y) as a function of position in the transverse plane. 
One sees that, for identical initial conditions, the HDM 
simulations lead to universally earlier freeze-out. This 
Figure also shows that for non-central Pb+Pb collisions 
(6 = 7 fm) the source at freeze-out is still larger in the y 
direction perpendicular to the reaction plane than in the 
reaction plane; the initial out-of-plane deformation of the 
reaction zone thus has not yet fully disappeared when the 
dynamical evolution was stopped. This explains why the 
elliptic flow had not yet saturated, see Fig. [l] and Fig. ^ 
below. 



C. Retuning the initial conditions for TTHM 

Our TTHM model does not include the hadronization 
of gluons and hence does not allow us to compute hadron 
spectra; hence we cannot directly compare the gluon mul- 
tiplicity density dN/dy and their /^-distribution from 
TTHM at £f to any data. However, the relation of 
these quantities to the observable hadron spectra after 
hadronization is not expected to depend on their dy- 
namical history prior to reaching the "freeze-out" point 
cf. Thus, if we want to create TTHM solutions which 
are likely to lead, after hadronization, to hadron spectra 
with similar normalization and shape as the data, we can 
use the gluon rapidity densities and fc^-spectra from ex- 
isting hydrodynamic calculations with initial conditions 
which were tuned to real data in , evaluate them at the 
same value of et where we stop the TTHM evolution, and 



± v ; 

FIG. 4: The gluon k± distribution from hydrodynamics 
(dashed line, with the parameters as given in the text) and 
TTHM (other lines). The TTHM curves correspond to 
the following sets of parameters: Dotted line: ro = 6fm/c, 
£o = eo/15; dot-dashed line: To = 7 fm/c, eo = eo/18; solid line: 
r = 7.8 fm/c, £ = e /20. Here e = 23.0 GeV/fm 3 is the ini- 
tial central energy density for 6 = Au+ Au collisions in the 
HDM. 



retune the TTHM initial conditions such that they repro- 
duce these HDM gluon spectra. With the TTHM spectra 
thus tied down to the final state in central (6 = fm) colli- 
sions from the HDM, we can then repeat our comparison 
of elliptic flow in the two models under more realistic 
boundary conditions. 

For the common decoupling point we chose in both 
TTHM and HDM the value e f = 0.35 GeV/fm 3 . This may 
appear somewhat low for still using only gluon degrees 
of freedom (with an ideal gas EOS in the HDM case), 
but it has the advantage that then the hydrodynamic 
model with the ideal massless EOS used here produces 
an elliptic flow V2 of roughly the same magnitude as that 
obtained for pions in Ref. Q using a more realistic equa- 
tion of state. Given the good agreement of those HDM 
calulations with the data Q, we can thus pretend to be 
comparing directly to data when comparing the TTHM 
results to our HDM reference. 

With this freeze-out energy density, our HDM refer- 
ence gives the gluon k± (or m±) spectrum shown in 
Fig. |4| which integrates to a gluon rapidity density at 
y = of dN/dy ~ 942. Given the large difference between 
the HDM and TTHM spectra for identical initial condi- 
tions shown in Fig. |2|, it is perhaps not surprising that we 
found it rather difficult to reproduce the HDM spectrum 
within the TTHM: one needs to find a way to dramati- 
cally lower the radial flow (since the freeze-out tempera- 
ture is fixed by Sf and is much higher in TTHM than in 
HDM) and at the same time increase the normalization 
of the spectrum in order to raise dN/dy. 

There are essentially only three parameters in TTHM 



FIG. 5: Variation of the gluon k± distribution with various 
TTHM model parameters. The solid and dashed lines repeat 
those from Fig. ^ for et= 0.6 GeV/fm 3 for comparison. The 
dot-dashed line is obtained by starting the system earlier at 
to — 0.1 fm/c. The long-dashed line has a lower initial energy 
density Si = eo /4. In each case the remaining parameters are 
held fixed at the values used in Fig. 0. 



that we can play with: eo, tq and £{. Examples of the 
effects of varying them are shown in Fig. |^. The lines 
from Fig. |^ corresponding to £f = 0.6 GeV/fm 3 are re- 
peated here for comparison. From Fig. || and Table | 
we learned already that changing et has no influence on 
dN/dy and little effect on dN/dyk±dk±. This leaves us 
only with changes of the initial conditions eo an d tq. The 
dot-dashed curve in Fig. || shows that decreasing t from 
0.6 to 0.1 fm/c goes in the wrong direction, by reduc- 
ing dN/dy without having much effect on the slope of 
the k± spectrum. Reducing the initial energy density eo 
by a factor of four produces the long-dashed line in the 
figure. Due to the lower initial energy density there is 
less time until freeze-out to produce radial flow, and the 
resulting spectrum is steeper, as desired; unfortunately, 
its normalization dN/dy decreases, too. Combining the 
insights from these two trial runs, we see that making 
the TTHM spectrum sufficiently steep requires still much 
smaller initial energy densities eo, combined with much 
larger starting times To for the TTHM transverse dynam- 
ical evolution, in order to increase dN/dy. This implies a 
long transverse thermalization time, which poses intrinsic 
consistency problems to be addressed later. 

Following this route, we find that the HDM gluon 
spectrum can be roughly reproduced with either one of 
the following three sets of parameters: To = 6 fm/c with 
£o = eo/15, To = 7fm/c and eo = eo/18, or To = 7.8fm/c 
and eo = eo/20. The corresponding spectra are shown 
in Fig. ^. These parameter sets are not unique and the 
agreement with the HDM spectra is not perfect, but what 
is clear is that a large reduction of eo by at least 1/15 
is essential in TTHM to obtain the much deeper HDM 
slope. 
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FIG. 6: V2 as a function of impact parameter from TTHM 
(solid) and HDM (dashed) simulations tuned to produce iden- 
tical angle- averaged spectra dN/dy k±dk±. The parameters 
used and notations here are the same as those in Fig. H. 
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FIG. 7: «2 as a function of fcx for two impact parameters 
(as indicated) from TTHM (solid) and HDM (dashed) simu- 
lations giving identical angle-averaged spectra dN/dy k±dk±. 



D. Elliptic flow from TTHM with retuned initial 
conditions 

Figure |^ again shows the impact parameter depen- 
dence of v 2 , but now with retuned TTHM initial con- 
ditions as described in the previous subsection. With 
the much lower eo forced upon us by the given slope of 
the final k± spectrum, there is much less time to generate 
transverse flow. Figure || shows that this leads not only 
to reduced radial flow, as reflected in the steeper single- 
particle spectrum for central collisions (see Fig. ||), but 
also cuts down on the elliptic flow v 2 which now remains 
significantly below the HDM level. As the latter is rep- 
resentative of the data [B O , we conclude that a model 
with only transversally thermalized momenta in the early 
collision stages cannot produce as much elliptic flow as 
required by experiment. 

In order to further strengthen this argument let us look 
at the k± dependence of v 2 and study (see Fig. [?]) 



/ d<P cos(2<£) 



v 2 (k±;b) = 



dN(b) 



(5.6) 



dy k± dk±d<p 



as well as its impact-parameter averaged ("minimum 
bias" §) value (see Fig. |) 



v 2 (k±) 



JbdbJ, 



dN(b) 



(5.7) 



/ k±dkj 



We here select for our comparison with HDM one of 
the above three sets of initial parameters for TTHM, 
namely tq = 7fm/c and Eq = eo/18, corresponding to the 
dot-dashed line in Fig. [|. 

Figures ^ and || show that the /c^-slope of v 2 is much 
lower than in the HDM. This is not only true for the 
impact parameter averaged elliptic flow shown in Fig. |^, 
but holds universally for all impact parameters; two ex- 
amples are plotted in Fig. [?]. For v 2 (k\^) from minimum 



bias events accurate data were published in which 
follow essentially the dashed HDM line in Fig. g| and def- 
initely exclude our TTHM. 

In Figure |{] we show the time evolution of the momen- 
tum anisotropy [^) 



(rpxx _ rpyy\^ 

(t ix + Tyy) 



(5.8) 



where the angular brackets indicate an average over 
the transverse plane. At freeze-out, this momentum 
anisotropy is translated into elliptic flow v 2 , where the 
coefficient between the two variables depends on the par- 
ticle rest mass ||. Figure ^| shows that with the hard 
equation of state for non-interacting massless gluons and 
RHIC-type initial conditions freeze-out actually happens 
before the momentum anisotropy saturates [ p"7| (and be- 
fore the initial spatial anisotropy has fully disappeared). 
This is even more true for the TTHM (with its even stiffcr 
EOS of P± = ±e) than for the HDM (where P = ±e): 
the TTHM freezes out significantly earlier, even though 
about the same amount of transverse flow is generated. 



E. The low-fcx limit of V2(k±) for massless particles 

A careful inspection of Figures ^ and |§| shows that at 
small transverse momenta the shape of the /c^-differential 
elliptic flow v 2 (k±) is different for TTHM and HDM. Fig- 
ure further shows that the small-fcx behaviour of the 
elliptic flow coefficient differs for massive and massless 
bosons. In Appendix [B] we present a detailed analytical 
treatment of the limiting behaviour of v 2 (k±) for small 
transverse momenta. For massive particles, the elliptic 
flow vanishes quadratically at k± — > 0, 



v 2 (k±) 



O(kl) 



(massive particles), 



(5.9) 




FIG. 8: Impact parameter averaged «2 as a function of k± 
from the two models. The dashed line is from HDM, the solid 
line from TTHM with initial conditions tuned to produce the 
same angle- averaged transverse momentum spectrum at the 
same value of ej. 



going over to an almost linear rise at higher k±_ (which 
eventually turns over to saturate at v 2 = 1 as fcj_ — > 00 
0). The transition from the quadratic rise at low k± 
to the quasi-linear behaviour at intermediate k± occurs 
around k± ~ m. In detail the transition between these 
regimes is complex and depends on the ratio m/T of the 
rest mass to the freeze-out temperature (see Fig. [l(]): for 
particles with m/T < 1 the elliptic approaches the curve 
for massless particles while the elliptic flow of heavier par- 
ticles with m/T 3> 1 remains always significantly below 
that of massless particles. 

For massless bosons, the singularity of the Bose distri- 
bution at zero momentum leads to a qualitative change 
in the low-momentum limit of the elliptic flow. For the 
TTHM, which lacks a thermal spread of the longitudinal 
momenta, the elliptic flow of massless bosons (gluons or 
photons) approaches a positive constant with finite pos- 
itive slope: 

v%(ki_) = a + bk±_ + 0(k 2 ± ) (TTHM, massless bosons). 

(5.10) 

In the HDM, the longitudinal thermal momentum smear- 
ing weakens the Bose singularity and causes v-i to still 
vanish at zero transverse momentum, albeit just barely 
so: as k± — > 0, V2(k±) vanishes with infinite slope: 

const 

V2(k±) = ~r~ (HDM, massless bosons). (5.11) 

m J— 

« J. 

It would be interesting to verify this prediction of the hy- 
drodynamic model for photons which are directly emitted 
from the expanding fireball. 



FIG. 9: The asymmetry measure e v as a function of the du- 
ration of the time evolution at b — 7 fm with similar final 
particle spectra from both models. Again the solid is from 
TTHM and the dashed line from HDM. 



VI. CONCLUSIONS 

We have tried to construct a macroscopic dynami- 
cal description of the fireball expansion in relativistic 
heavy-ion collisions which relaxes the usual assumption 
of complete local thermal equilibrium which underlies 
the popular and successful hydrodynamical approach. 
This addresses, on a macroscopic level, the observation 
that existing parton-based microscopic kinetic models 
p(j| . [2]] , p2| seem unable to generate sufficiently fast ther- 
malization to reach the hydrodynamic limit. We have 
tested a toy model where transverse momenta are ther- 
malized very quickly while the system is streaming freely 
in the longitudinal direction (i.e. without longitudinal 
momentum transfer). The initial longitudinal motion 
was assumed to be boost-invariant, and the model pre- 
serves this property dynamically at later times. 

Starting from the kinetic transport equation, we de- 
rived the corresponding macroscopic equations of motion 
(TTHM) and solved them in parallel with the usual hy- 
drodynamic (HDM) equations. Since our TTHM model 
is parton-based and lacks a description of hadronization, 
we decided to compare the two approaches just before 
the onset of hadronization. To make phenomenologically 
meaningful comparisons we used a HDM reference which 
had previously been shown to yield a good description of 
RHIC data when allowed to evolve until hadron freeze- 
out. 

The TTHM turns out to have an extremely stiff equa- 
tion of state, P± — |e, which causes severe phenomeno- 
logical problems: in order to avoid the creation of too 
much radial transverse flow, which would render the 
single-particle spectra much flatter than observed, while 
preserving the total multiplicity, we had to start the 
TTHM evolution very late, at already quite low energy 
densities. This implies a very long transverse thermaliza- 
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FIG. 10: V2 as a function of k±_ for massless and massive 
bosons from the two models. The solid lines are for massless 
bosons while the dashed and dot-dashed lines are for bosons 
with m=100MeV and m=l GeV, respectively. 



tion time scale (and, of course, an even longer longitudi- 
nal one) . While this cannot a priori be excluded theoret- 
ically, it causes consistency problems with our implemen- 
tation which assumes that transverse collective dynamics 
sets in only after we start the TTHM evolution (i.e. only 
after thermalization). This is not realistic: even before 
complete thermalization of the transverse momenta the 
system will start to evolve in the transverse directions, 
since nothing confines it. This pre-equilibrium transverse 
dynamics will likely have a collective component affecting 
the transverse momentum spectra, necessitating an even 
further shift of the starting point for TTHM evolution. 

We have not followed this issue any further, since there 
is also a severe problem when comparing the TTHM pre- 
dictions with the elliptic flow data, which in our opinion 
eliminates it as a viable model. For the same radial flow, 
constrained by the slope of the single-particle spectra, 
TTHM provides only about half as much elliptic flow 
as the hydrodynamical model and thus significantly less 
than that required by the RHIC data g g, g, §. This is 
the central result from the present study, and it strength- 
ens the previously made argument || [IT] that the large 



elliptic flow measured at RHIC can only be reproduced 
by models which assume very rapid thermalization of the 
momentum spectra in all three dimensions. Partial ther- 
malization, such as the effectively only two-dimensional 
thermalization in the TTHM model, is not enough to gen- 
erate the observed elliptic flow. This leaves us with the 
continued challenge of identifying the microscopic mech- 
anisms within QCD which can cause such rapid three- 
dimensional thermalization. 
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APPENDIX A: DECOMPOSING THE 
ENERGY-MOMENTUM TENSOR 

To calculate the coefficients in the decomposition 
( |2.19| ) of the energy-momentum tensor we find it easi- 
est to use th e form ( 2.12 ) for the distribution function 
in Eq. ( |2.1g| ). Assuming massless particles (gluons), we 
write the four-vector fc M as 



fc M = k± (cosh y, cos 0, sin 0, sinh y) . 



(Al) 



We integrate over d 3 k/ko — dy k±dk± d<j) in Eq. ( 2.18 ) by 
first doing the trivial 77-integration and then integrating 
over k±. After int egrat ion over 77 the exponent of the 
Bose factor in Eq. (2.12) reduces to 



k-u l±k±, f . 7xfcx,. , 
— = — jr- (1 - k±-v±) = — y~ C 1 - v ± COS{(j)~(p s )) 

(A2) 

where <p—(p s is the angle between k± (azimuthal angle </>) 
and v± (azimuthal angle <f> s ). The fcj_-integration is then 
easily done using 



k n dk. 



n\((n+ 1) 



- 1 



(A3) 



In our case a = 7_i_(l — v± cos(0— <j) s ))/T and n — 3, and 
the relevant zeta function is £(4) = 7r 4 /90. In this way 
we arrive at 



T V TT 2 T 4 
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(A4) 



12 



/ cosh 2 77 cos (j) cosh 77 sin </) cosh 77 cosh 77 sinh 77 \ 



cos <p sm i 

,2 / 



cos 4> cosh 77 cos (/> 

sin cosh 7/ cos <fi sin sin" 
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cos (/> sinh 7/ 
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sinh 2 77 / 



(A5) 



The remaining azimuthal integration is done by shifting the integration variable to 9 = < 



using 



in the numerator M F , and exploiting the formula 

r7r dO 1 1 



;Pn- 



_ n 2ir (l-v± cos 6<) n (1-t; 2 )"/ 2 
where the P n are Legendre polynomials. We need 

7x-Pi(7x)=7x, 7i^ 2 (7±)-7i(l + M) 



= 72^-i(7x), 



7l^3(7x)=7i(l+l«i) 



M , M = M , and M have no 0-dependence and the corresponding angular integrals thus give 
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l_ n 2ir(l-v x cos(0-0 s )) 4 
For M 01 = M 10 and M 13 = M 31 we need 
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(A10) 



In the last step on the first line we dropped the vanishing term ~ sin 9 and also added and subtracted a term l/v± in 
the numerator. In the second line we used cos(f> s =v x /vj_. For the remaining components we proceed similarly. The 
energy-momentum tensor thus takes the form 



TO 7T 2 T 4 

r 9 60 



7i 



/ (1 + |?j 2 ) cosh 2 77 ^v x (4 + t; 2 ) cosh 77 \v y (A + 7jj_) COSI177 (1 + |«^) cosh 77 sinh rA 
±7^(4 + ?j 2 J cosh?; \{1 + Av 2 x - v 2 ) %v x v y 



hv y {A + v 2 ^) cosh 77 



\{l-v 2 x +Av 2 ) ± Vy (A + vl) sinh 77 



\(1 + \v\) cosh 77 sinh 77 \v x {A + t; 2 ) cosh 77 ^v y (4: + v 2 ^) sinhi] (1 + |7; 2 ) sinlr 77 / 



Contracting this result once with each of the vectors 7i u , n", and m u we obtain 
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A second contraction then yields the coefficients A through H in Eq. ( 2.19 ) as given in equations ( 2.22 )-( 2.27 ). 



APPENDIX B: THE SMALL TRANSVERSE 
MOMENTUM LIMIT OF ELLIPTIC FLOW 



In this Appendix we derive the small-fc^ limit of the 
elliptic flow coefficient v-i for TTHM and HDM, for mas- 
sive and massless particles. Suppressing the impact pa- 
rameter dependence which is irrelevant here, we write 



Eq. (5J3) in the form 



V 



(Bl) 



where the numerator TV and den omi nator T> are obtained 
from the Cooper- Frye formula (5.1). For massive parti- 
cles we must generalize Eq. (5.3) by writing the first term 
as mi cosh(?7— 77) and use 



k ■ u = 7_l (rn± cosh(?;— 77) — k±-v±) . 



(B2) 
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For the T THM the distribution function is given by 
Eq. (2.12), and we obtain for particles at midrapidity 
y = 



2. Massless bosons in the TTHM: For massless bosons 
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where 



ckf) ( cos(2</>) 
~2~i 1 1 
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g~f±(m± cosh7i — k±-v±)/T ^ ' 

(B4) 

Here 7j_ and Tf are functions of the position in the trans- 
verse plane, 7_l — j ±(x± ) and T{ = Tf(x±). For the HDM 
we instead use Eq. (3.1) and get 



d 2 x\ 



Tf 



T><h 



(B5) 



The main difference (except for the different weighting 
of the integration over the transverse plane d x± = 
r dr d<f> s ) is the nontrivial 77-integration in the HDM case. 

1. Massive particles: For non-zero rest mass, m ^ 0, we 
can expand in both cases the </>-integrand for small k±. 
Keeping only terms up to first order in k± we find 



m cosh 77 



T>di I ~ e -y±mcoshr]/T _ ^ 

k±-V±Tf\ ( ^ kj_-v±_ 



2 , ri <«> 



i - 



m cosh rj 



1 



7i 



^J 1 2 g — 7_L m cosh 7?/T 



Higher orders of k± come with higher inverse powers of 
T or m cosh 77; the corresponding ^-integrals thus are all 
finite (trivially so in the TTHM), and the k± -expansion 
is well-behaved. 

As in Appendix A we write kj_-vj_ = kj_v± cos(0— <j) s ) 
= fcj_t>j_ cos 9 and use — n < 9 < n as integration variable. 
Then, using Eq. @, 



k± ■ Vj_rf = (d r Tf) cos 8 - 



bTf ) sin < 



(B7) 



where d r Tf — v±-\7±Tf and d^Tf = z-(v±x\7 ±Tf). Here 2; 
is the unit vector along the beam direction, and v±_ co- 
incides with the radial unit vector, v±_ = (cos cj> s , sin <p s ). 
We also need 



cos(20) = cos(2</> s )(cos 2 



sm 



2 sin(2<^ s ) sin 8 cos 9. 



(B8) 



After th ese manipulations the angular integrations in 
Eq. (B6) are easily performed. Due to symmetric in- 
tegration over the full circle, only terms containing even 
powers of cos0 and/or sin0 survive. We find 

Nj, = 0{k 2 _), V 4> = a + bk ± +(D{kl), (B9) 

where a and b are integrable functions of x± and r\. As 
a result, for massive particles V2(k±) = 0(k 2 ), i.e. the 
elliptic flow coefficient vanishes with zero slope when 
k± — ►(), as first observed by Danielewicz fl9t]. 



at y = the angular integrals Af$, D§ in Eq. (B4) reduce 
to 



t_ r # r cos(2^) I i-^Str 

7-L 2vr \ 1 J 1 _ 

cosh 77 

(cosh 77 - k±-v±) 
X : . 

g iL_J~(coshJ7— fe-L-fx) _ ^ 



(BIO) 



We can try to expand the last factor for small k± by 
using 



^T = l-| + ^ 2 ). 



(Bll) 



However, since in this expression the expansion parame- 
ter is x = 7 ^" J - (cosh 77— k±-v±), each additional power of 
k±_ brings in another factor cosh 77. This is no problem for 
the TTHM with its factor 15(77) under the 77-integral, but 
for the HDM the 77-integrals of the expansion coefficients 
diverge, rendering the expansion meaningless. For the 
HDM thus the 77-integration must be performed before 
we expand for small k± (see 3. below). 

For the TTHM, on the other hand, we can continue 
along this direction. Using the same manipulations as in 
1. above, we find up to linear order in k± 



T r dd f cos(20 s ) cos(26>) - sin(20 5 .) sin(20) 1 

~]-^\ 1 J 

1 — (d r Tf) cos 9 — (d^Tf) sin 9 



x 1 



1 — 7Jj_ COS ( 
7-L&J 



2T 



-(1— v± cos( 



(B12) 



We use Eq. ( A7) to derive the following table of integrals: 
<•* d6 1 



.„. 27T 1 — v±_ COS 9 

w d8 cos6> 
_n 2ir 1— v±_ cos 9 
71 d8 cos 2 8 

_ w 27T 1 — v± COS 9 

* d9 sin 2 6» 

_ w 27T 1 — v± COS 9 

* d9 cos(26>) 
. w 2tt 1— v± cos 9 

71 d9 sin 2 9 cos 9 _ 
_ v 2tt 1—v± cos 9 

77 d9 cos(26») cos 9 
2n l—v± cos9 



IX, 
7-L-l 

1 

7J_l 

< ' 

7-L-l 

2 ' 

7-LUl 

1 

7± 
1 

2^T 



7-L-l 
v± 

7_l-1 



7±wj 
1 fl±~^ 



J±v± V v± 



(B13) 
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and find 



T 



7-1- 



cos(20 s ) I 1 - 



•sin(20 s ) 



d</>Tf 
l±v± 



T 1 



5 r r f - — 

7_lW_l / 2 



0(k] 



-O(kl). (B14) 



Now both the numerator and denominator approach 
nonzero constants as — > 0. Since the fireball center 
freezes out later than its edge, the term ~ (d r T(), after 
integration over the transverse plane, contributes with a 
negative sign, and the fcj_-independent first terms in TV^ 
and 2?0 are correspondingly positive. As a result, the 
small-fcx expansion of the elliptic flow coefficient takes 
the form 



V 



1 -bk± 
abk± + 0(kl), 



(B15) 

where a and 6 are positive constants arising from the 
integration over the transverse profiles of the freeze-out 
time and flow velocity. Thus, as k±_ — > 0, t>2 approaches 
a positive constant value with nonzero positive slope, as 
confirmed by Figs. @ to [)[ 

3. Massless bosons from the HDM: In the last step we dis- 
cuss the corresponding limit for massless bosons from the 
HDM, in order to understand the surprising behaviour 
shown in the upper part of Fig. ^[ As already mentioned, 
in this case the ry-integration must be performed before 
expanding in k±. Doing this integration first, we rewrite 
Eq. (PI as 



r dr d(f>i 



Tf(r, 4> s 



(2ir) 2 J ~" r ° 7x(r, <fr s 
d6 J cos(2<£ s ) cos(26l) - sm(2(/) s 
2^ 1 1 



sin(20) 



dr] 



£ ^coshr; 



k±-V±Tf 



where all the dependence on the magnitude of fcj 
den in the variable 



T 



(B16) 



is hid- 



(B17) 



Let us denote the result of the 77-integral by I v . Expand- 
ing the Bose-Einstein distribution in a power series and 
performing the ry-integration term by term we find 



In = 



r ( >•>• i '■■■>'■ I dr}e~ n ^"" u '' 



x (cosh ry — (d r T{) cos 9 



jTf) sin ( 



2^, 

71=1 



,nv±£ cos 9 



■(cos6>(9 r r f ) + sm6/(^r f ))Ko(nO , (B18) 



where the K„ are modified Bessel functions of the second 
kind. The integration over the momentum-space angle 
9 can now be performed, too, yielding modified Bessel 
functions of the first kind I„: 



r d9 

/ COS! Vb 

J-* 27T 



(B19) 



where z = nv±£. After some simple manipulations us- 
ing trigonometric identities we find for the elliptic flow 
coefficient 



v 2 (k±) 



EZJrdrd^^)V n (r,, 



(B20) 



7±{r,<t>s ) ' 
Tt(r,(f> 3 

where 

W n = cos(20 s )[ei 2 (n^OKiK) 

-^Y e(li(n»xO + Is(n»xO)Ko(ne)" 

+ sin(20 s ) ^ i (l x (nv ± - I 3 (n«x£)) K o«) - 

T> n = CIo(nwxOKiK)-(9 r T f )CIi(nwxOKoK). 

(B21) 

These expressions generalize the result given in [ ?J to the 
general case of a freeze-out time Tf which depends on the 
transverse position x±. 

We can now study the small-momentum limit of v 2 by 
letting £ — > 0. Let us first consider the Boltzmann ap- 
proximation which corr esponds to keeping only the terms 
M\ and T>\ in Eq. (B20). Using the small-£ expansions of 
the Bessel functions one finds that in this approximation 



lim v- 



Boltz 



In — 

T 2 T 



(B22) 



For massless Boltzmann particles, the elliptic flow from 
the HDM thus goes to zero with zero slope, but by a 
factor hi(k±/T) more slowly than for massive particles, 
see point 1. above. 



For bosons we must sum over all n in Eq. ( |B20|) . For 
k± — > 0, i.e. £ — > 0, the sum over n becomes the discrete 
representation of an integral: defining 



we have 



/^(C) = l>xC) M0> 



(B23) 



n=l 



€/2 



/»x(C)dC- 



(B24) 



By expanding the Bessel functions for large arguments 
one sees that for large £ 



1 



-(l-«-L)C 



(B25) 
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such that the integral in Eq. (B24) clearly converges at 
the upper end. The fc^-dependence of V2 can now be 
extracted by studying the dependence of this integral on 
its lower limit £/2 = j±k±/2T. To do so we employ 
the following trick: We split the integral at an arbitrary 
intermediate scale A, 



fvAOdC = / fvA()d(+ / fvAOdC 
e/2 Jap J\ 

A 

f v± (()dC + F(\;v ± ), (B26) 

e/2 



where F(X;v±) is a finite number whose dependence on 
A cancels against the A-dependence of the first integral, 
irrespective of how we choose A. By taking £/2 < A <C 1 
we can evaluate the first integral analytically by making 
use of the expansion of the Bessel functions for small 
arguments, keeping only the leading terms. In the limit 
£ — ► we find in this way 



£/2 



Ia(vi0Ki(O dC 
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«/2 
A 

A 

€/2 



l!(i;xOKo(C)dC 



^ln«) + C7(A), 



I 3 KC)K (C) dC^O- 



ioKC)Ki(C)rfC 



hxj+C(X), (B27) 



where a is a A-independcnt constant while C, C de- 
pend on A but combine with the corresponding constants 
F(X;v±) to finite, A-i ndep endent con stants. Inserting 
these results into Eqs. ( B2C ) and ( B21 ) we finally obtain 



const. 

hm V2[k±) = - — Tfr- 



(B28) 



which approaches zero with infinite slope. This is ex- 
actly the behaviour seen in the upper panel of Fig. ^ for 
massless bosons. 
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